NMR and LDA evidence for spiral magnetic order in the chain cuprate LiCu 2 02 
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We report on 6,7 Li nuclear magnetic resonance measurements of the spin-chain compound LiCu2 02 
in the paramagnetic and magnetically ordered states. Below T~24 K the NMR lineshape presents a 
clear signature of incommensurate (IC) static modulation of the local magnetic field consistent with 
an IC spiral modulation of the magnetic moments. 7 Li NMR reveals strong phason-like dynamical 
fluctuations extending well below 24 K. We hypothesize that a series of phase transitions at 24.2, 
22.5, and 9 K reflects a "Devil's staircase" type behavior generic for IC systems. LDA based 
calculations of exchange integrals reveal a large in-chain frustration leading to a magnetical spiral. 



Despite extensive efforts during the last decades spin or- 
dering in frustrated 5=1/2 quantum spin chains still 
remains a matter of broad activities. Q, [3- EJ Dl 0- El 
Rich phase diagrams with commensurate (C) and incom- 
mensurate (IC) phases, with spin- and charge ordering, 
dimcrization, or superconductivity have been predicted. 
Most studies have been focused on various cuprates with 
corner- or edge-shared Cu04 plaqucttes. Edge-sharing of 
CuC>4 plaquettes leads to Cu02 chains with a nearly 90° 
Cu-O-Cu bond angle causing a reduced nearest neighbor 
(nn) transfer and a next-nearest neighbor (nnn) trans- 
fer of similar size allowing frustration effects. IC spiral 
states driven by ferromagnetic (FM) nn exchange and 
in-chain frustration have been predicted theoretically for 
Cu02 chain compounds such as Ca2Y2Cu50io but dis- 
carded experimentally. 0, 0| The observation of in-chain 
IC effects in undoped quasi-lD cuprates is so far re- 
stricted by a sharp magnetic field driven C-IC transition 
observed in the spin-Peierls system CuGe03 at high mag- 
netic fields. Here we report on 6,7 Li nuclear magnetic 
resonance measurements (NMR) and local density (LDA) 
based analysis of the electronic and magnetic structure 
of the chain compound LiCu202- We show that the ob- 
served spontaneous magnetic order can be described by 
a spiral modulation of the magnetic moments. Indepen- 
dently, LDA calculations and a subsequent Heisenberg- 
analysis reveal strong in-chain frustration driving spiral 
ordering in accord with the NMR data. 



LiCu202 is an insulating orthorhombic compound[_ 
with bilayers of edge-shared Cu 2+ -0 chains running 
parallel to the 6-axis separated by Cu 1+ planes. It ex- 
hibits a high-temperature antifcrromagnctic(AFM)-likc 
Curie- Weiss susceptibility v(T). Low-temperature x(^) 
and specific heat stufies [lj, [llj point to a series of in- 
trinsic phase transitions at Tw24.2 K, T«22.5 K, and 
T«9 K pointing to a complex multi-stage rearrangement 
of the spin structure. Magnetization studies performed 



in external fields up to 5 T did not reveal any signa- 
tures of field-induced transitions. /iSK data[ll| point 
to a broad distribution of magnetic fields at the muon 
stopping sites. LSDA(+U) calculations point to an FM 
in-chain ordering. However, a simple FM ordering is 
in conflict with the fj.SK data[llj and the AFM dimer 
liquid picture. P, Thus, to settle origin and charac- 
ter of the spin ordering in LiCu202 further studies are 
required. NQR and NMR measurements being a local 
probe are standard tools to elucidate the electronic and 
the magnetic structure. However, the Cu NMR signal in 
LiCu202 stems presumably from nonmagnetic Cu 1+ ions 
[l3| and reflects mainly the disorder due to Li-Cu 11011- 
stoichiometry to be discussed elsewhere. We have per- 
formed 6,7 Li NMR measurements on large single crystals 
at several temperatures both in the paramagnetic [l3j 
and in the ordered phases. The spectra were measured 
for two main orientations: H || (a, b) and H _L (a, b), 
and by 10° steps inbetween, by sweeping the magnetic 
field at a fixed frequency of 33.2 MHz, the signal was 
obtained by integrating the spin-echo envelope. In the 
paramagnetic state (T=45K) the 7 Li NMR line (Fig. [TJ 
has a full width of 0.01 T (at the basement) and shows a 
typical (for a 7=3/2 nucleus) first order quadrupole split- 
ting pattern with vq=51.7 kHz. For T<24 K a dramatic 
change of the 7 Li NMR line shape is observed. In general, 
the low-T NMR line shape shows a wide continuous sym- 
metric distribution which evolution is presented in Fig. 
121 for the H || (a, b) geometry. In the ordered phase we 
still deal with a strong and narrow central peak but with 
a well developed inhomogeneously broadened basement. 
Lowering T the central peak weakens and disappears near 
TwlO K, while the basement gradually evolves into a well 
defined sextet whose full width exceeds by a factor of 30 
(!) that of the paramagnetic phase. A similar situation 
takes place in H 1 (a, b) geometry, but, with a spe- 
cific quartet shape of the low-T saturated NMR spectrum 
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FIG. 1: The 7 Li NMR spectrum at T=45 K for different 
orientations of the magnetic field with respect to the c-axis. 
The solid curve on the upper spectrum: our simulation. 



(Fig. |3J) . Such spectra are unique signatures of an infinite 
number of non-equivalent Li sites. They correspond to an 
IC static modulation of local fields |l5lll6U17l | To discrim- 
inate between magnetic and electric (quadrupole) effects 
we have performed parallel 6 Li NMR studies. |l8| The 
puzzling closeness of two NMR spectra below the critical 
temperature (Fig. |3| top panel) unambiguously points to 
a distribution of local magnetic fields as the origin of the 
observed phenomena. The additional local magnetic field 
h(R) seen by a nucleus at a site R near a CuC>2 chain can 
be directly converted to a local spin polarization on the 
neighboring sites S(R + r): h(R) = £ r ^(r) s ( R + r), 
where A(r) is the anisotropic tensor taking account of 
magnetic dipole and (super)transferred Li-Cu(O) hyper- 
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FIG. 2: The 7 Li NMR spectra taken at various temperatures 
below the critical temperature of the magnetic ordering. The 
magnetic field is directed || to the a(b) axis. 
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FIG. 3: The 6 Li and 7 Li NMR spectra taken at 4 K. The 
solid curves in both panels show the results of our theoretical 
simulation (see text). 



fine interactions. When h <C Hq (the uniform external 
magnetic field: Hq = i?o(sin cos $, sin O sin<!>, cos 0)): 
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It should be noted that the IC struc 



4>2 are phase shifts, 
ture can be described in different ways, both as a simple 
plane-wave or as a regular array of domain walls (soli- 
tons). If the concentration of solitons is large enough, 
we arrive at a soliton lattice. Below we assume a plane- 
wave modulated IC spin state, where the spin modulation 
phase is a harmonic function of the space coordinate in 
the direction of the modulation (hereafter Z-axis). Such 
a spiral spin structure we define as follows: S(R) = 
S^sin 0(R) cos </>(R) , sin 0(R) sin 0(R), cos 0(R)), where 
0(R) = fciZ+0o, <MR) = ^2^+00 with fci, 2 being IC spa- 
tial modulation wave numbers. In particular cases this 
spin structure can be reduced to any types of C states 
including simple FM- and AFM ones. In the continuum 
approximation the resultant NMR line shape for a static 
case corresponds to the density of the local field distri- 
bution g(Z) oc \dH]Z)/dZ\~ 1 convoluted with the indi- 
vidual line shape. Actually the straightforward ap- 
plication of this procedure is restricted by rather simple 
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particular cases. However, the lack of a central peak 
and of any other sharp lines in the NMR spectrum at 
T=4 K for all orientations of the external magnetic field 
points to a spatial modulation of all three spin com- 
ponents. Hence, we are led to consider complex spa- 
tial spirals with a modulation of the longitudinal (Sz) 
as well as of the transversal (Sx,y) spin components. 
The specific quartet line shape for the H _L (a, b) ge- 
ometry and the sextet one for H || (a, b), with equiv- 
alent signals for Ho || b, Ho || a, and the field ori- 
ented along the diagonal in afc-planc, yield constraints 
for possible types of spirals reducing substantially their 
number. The measured NMR spectrum can be well de- 
scribed by a simulation assuming a spin spiral twisting 
along the chain b axis with identical modulation periods 
both for the longitudinal and the transversal spin compo- 
nents (ki=k2). Within the long-wavelength limit and for 
a short-range hyperfine coupling the NMR field shift due 
to a (super)transferred hyperfine coupling in H 1 (a, b) 
geometry for the spiral with a "positive" twisting along 
b-axis (k || b, or k > 0) can be written as follows: 
h = (5/4) (E^xx) [sin(2/cZ+0i ) -sin <f>i] , where S £ A xx 
is the maximal value of the hyperfine field for the collincar 
spin ordering assumed. The spectral density function 
g(H) has a rather simple form g(H) oc [1 — (w(H) + 
sin0i) 2 ]~5, with w(H) = 4{H - H )/SA XX . As a result 
we arrive at two singularities uj±(H) = ±1 — sin^i with 
a separation \H + — H- \ = (S/2)A XX in between, shifted 
with regard to Hq by {S/A)A XX sin0i. Taking into ac- 
count both "positively" (k > 0) and "negatively" (k < 0) 
twisted spirals we arrive at a specific symmetric quartet 
NMR line shape. The situation looks slightly more com- 
plex for the "in-plane" geometry, where we deal with two 
types of twins. For Ho || b or a axis we arrive at a su- 
perposition of two different NMR signals stemming from 
differently oriented twins, with Ho || Z or Ho _L Z, re- 
spectively. The former are described by a simple expres- 
sion for the NMR frequency shift: h = (S/2)A ZZ cos kZ 
with the well-known symmetric doublet line shape [lij . 
while the latter yield a NMR response qualitatively sim- 
ilar to that of H _L (a, b) geometry. Thus we arrive at 
a sextet line shape. Making use of the Gauss convolu- 
tion procedure we successfully simulate the NMR 
spectra both for H _L (a, b) and H || (a, b) geometry 
(Fig|3J). It is noteworthy that the current NMR study of 
IC phase does not yield information regarding both the 
IC period and the absolute value of the spin order pa- 
rameter. Nevertheless, the simulation provides valuable 
insight into the effective (super)transferred Cu(0)- 6,7 Li 
hyperfine interactions. In fact, we obtain: SI^A^I ss 
S\J2^x\ «3 kOe, S\Y,A y y\ «2.6 kOe. The large 
magnitude and closeness of these parameters point to 
the predominance of covalently enhanced isotropic (su- 
per)transferred Li-Cu(O) hyperfine interactions as com- 
pared with anisotropic magnetic dipole terms. It should 
be noted that our current approach implies implicitly a 



mutual collinear interchain ordering. A more detailed 
analysis of interchain coupling and impurity effects are 
challenging problem left for future studies. 

The discreteness of the crystal lattice causes pinning 
which tends to lock in the IC wave vector for a spiral at C 
values giving rise to a "Devil's staircase" type behavior. 

We hypothesize that LiCu2C>2 with a series of phase 
transitions at 24.2, 22.5, and 9 K exhibits such a be- 
havior. This conjecture is supported by the temperature 
evolution of our NMR spectra (see FigEJ which reveals a 
well developed central peak in a wide temperature range 
below 24 K. This is a signature of strong low-energy dy- 
namical fluctuations of local magnetic field typical for 
acoustic-like phase fluctuations, or phasons generic for 
IC structures. So, the first order phase transition at 9 
K may be attributed to pinning effects. Noteworthy the 
phason-like spin fluctuations may result in a sizeable fre- 
quency dependence of the transverse relaxation rate T^" 1 
affecting strongly the NMR line shape. These effects dif- 
fer substantially for different singularities in the NMR 
line shape, in particular, they maybe responsible for the 
weak intensity of the outer doublet in the NMR spectrum 
for H || (a, b) geometry. 

To elucidate the origin of the observed IC spin order- 
ing, the knowledge of the relevant magnetic couplings 
is essential, especially with respect to frustrations. In 
particular, the ratio a=\J nnn / J nn \ in a single-chain ap- 
proach of the FM-AFM J\- J2-Hciscnberg model is deci- 
sive for C order (a < 1/4) or IC one (a > 1/4). To esti- 
mate the main transfer integrals t and exchange couplings 
J, we performed full potential LDA calculations using the 
FPLO code.^^ The resulting bands near the Fermi en- 
ergy, shown in Fig. 0] have been fitted by an extended 
tight-binding (TB) model to obtain the Vs. The result 
is given in Tab. [I] The in-chain dispersion along the k y 
directions (r-Y) shows deep minima in the middle of BZ 
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FIG. 4: Band structure of L1CU2O2 near the Fermi energy, 
with r=(0,0,0), X=(1,0,0), Y=(0,1,0), Z=(0,0,1), R=(l,l,l), 
S=(1,1,0), T=(0,1,1), U=(1,0,1) in units of 7r/a,7r/6,7r/c, re- 
spectively. The inset shows a sketch of a bilayer with double 
chains (• and filled O) running along b together with our 
notation for the main transfer and exchange processes. 
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TABLE I: LDA derived transfer integrals t and exchange 
constants J in meV, notation according to inset of Fig. [I] 



indicating the predominance of the nnn transfer over the 
nn term with i2y~-109 meV and \t y \~ 64 meV, respec- 
tively. This causes a strong in-chain frustration within 
the spin- 1/2 Heisenberg model H = -A .S.S ,. The 

AFM contributions to the J's estimated within a single 
band Hubbard model are given by J^ FM = At^/U. We 
adopt for the screened on-site repulsion U=3 eV as an 
appropriate value ||. The result is given in Tab.[I] The 
screened FM Heisenberg contributions J FM (see Tab. [|J 
have been found using Wannier functions [2(J of Li 2 Cu0 2 
which is closely related to LiCu202 with respect to its in- 
chain geometry and crystal field. The effective exchange 
j e ff i s given by the sum of the two former contributions 
(see Tab.[|J. As one may expect for a nearly 90° Cu-O-Cu 
bond angle according to the Goodenough-Kanamori rule, 
we arrive at a neg ative (FM) total nn in-chain exchange 
Jj,«— 8.1 meV,|2lj and at an AFM nnn exchange: J-2 y ~ 
+14.5 meV. The main interchain transfer process t x is 
comparable with the nn term t y , but the corresponding 
FM exchange contribution is small resulting in J**> =5.7 
meV. For the small transfer terms t xy and t x , the corre- 
sponding FM contributions are expected to be very small 
and are ignored therefore. From the values in Tab.[l]wc 
conclude that LiCu2C>2 is a quasi-2D system in an elec- 
tronic sense but magnetically more anisotropic. Using 
the calculated J's, we end up with ct~1.8 well above the 
critical ratio of 1/4. Wc emphasize that for a 3d ar- 
rangment of chains such as in LiCu202 and the FM sign 
of J y found, the in-chain frustration is the only strong 
enough source to drive a spiral. Thus LDA provides 
strong microscopic support for the scenario suggested by 
the phenomenological NMR analysis. In a more realistic 
anisotropic multiple chain situation, a c would be slightly 
larger. However, with a~1.8 the ignored interchain cou- 
pling and spin anisotropy should not change essentially 
our physical picture. This is at variance with the ho- 
mogeneous FM in-chain order found in Ca 2 Y 2 Cu50io 
(a = 0.05) and Li 2 Cu02 (a< 1) with much stronger 
interchain coupling. Hence, LiCu 2 02 is unique, in be- 
ing the first realization of a long sought for spiral chain 
compound. It is the most strongly frustrated cuprate we 
know. For such a compound the study of doping and in 
particular the search of unconventional superconductiv- 
ity are intriguing issues of general interest. 

In conclusion, we present the first experimental evi- 
dence of a low-temperature IC in-chain spin structure in 



an undoped quasi-lD cuprate. |22j Below 24 K the 7 Li 
NMR lincshape in LiCu202 presents a clear signature 
of an IC static modulation of the local magnetic field 
consistent with a spiral modulation of the magnetic mo- 
ments twisted along the chain axis. The intrinsic incom- 
mensurability accompanied by pinning effects is believed 
to be responsible for an unconventional evolution of a 
low-temperature spin structure with a series of succes- 
sive phase transitions resembling the "Devil's staircase" 
type behavior. Applying an independent extended TB 
and subsequent Heisenberg analysis based on LDA calcu- 
lations we succeeded to elucidate a strong in-chain frus- 
tration and this way also the microscopic origin of the 
magnetic spiral seen by the NMR. 
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